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Abstract 

We perform a model study on the 2-flavor lattice Schwinger model 
using standard Wilson and 0{a) improved Sheikholeslami-Wohlert 
(SW) action. We find, that the phase diagram is altered, the crit- 
ical line shifted closer towards its continuum value Kc = 0.25. We 
find no improvement in the rotation invariance of meson propagators; 
the scaling of the Schwinger mass is considerably improved, high mo- 
mentum states are not. The additional cost of « 30% CPU-time is 
highly justified when calculating masses, but not for high momentum 
observables. 
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1 Introduction 



In recent years, there has been an increasing interest in possible use of im- 
proved actions to lattice QCD. There exist different improvement schemes for 
lattice actions and observables, all which aim at the extraction of continuum 
physics from simulations on coarse lattices. 

The simplest of these improvement schemes is the clover action, an 0{a) 
improvement of the Wilson fermion action suggested by Sheikholeslami and 
Wohlert who implemented the Symanzik improvement [2 for the Dirac 
operator. Several calculations in quenched and some in full QCD have been 
performed using this action especially to calculate both heavy and light 
hadron masses. The results obtained seem promising, so that one might 
want to check the improvement effects of the clover action in a simple and 
well-known toy model. Here we study the effects of the 0{a) improvement 
of the Wilson action in the full lattice Schwinger model. This model seems 
optimally suited for this purpose. 

Massless QED in 1 + 1 dimension was first studied by Schwinger [H] as an 
example of an explicitly solvable QFT. In subsequent papers |3] the model 
was generalized to allow for fermion masses and include different flavors. 

The Schwinger model became of considerable interest during the early 
1970's, because it shows some remarkable features - like fermion confinement 
or a G- vacuum structure - which are reminiscent of QCD. For the purpose 
of this paper, it is sufficient to note, that the massless A^-fiavor Schwinger 
model has a spectrum which consists of one massive isosinglet boson and a 
isomultiplet of (A^^ — 1) massless fermion-antifermion bound states (for fur- 
ther discussion cf. e.g. |^). The mass of the isosinglet state, the 'Schwinger 
mass' (in lattice units) is given by 



where /? = l/(e^a^) denotes the dimensionless gauge coupling. 

In the next section some technical details about the definition and the 
numerical simulation will be given and in Sec. 3 we discuss the effect of the 
improvement on the rotation invariance of the propagator, the phase diagram 
and the dispersion relations of the Schwinger boson. 
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2 0{a) improvement and simulation 



The fermionic part of the Wilson action for the massive, N-flavor Schwinger 
model is 

X X,fl 

The fields have 2N components, o"^ are the Pauli matrices. This action has 
a classical 0{a) error introduced by the Wilson term, which can be removed 
by rotating the fermion fields 

^ (1 _ 1^)V1/ , ^r^fr(l + i^) (3) 

and adding a term 

1 _ 

Ssw = -i^cswt: (^) 

^ X 

to the action. Ff^^, denotes the clover-leaf discretization of the field strength 
tensor. The resulting lattice Dirac operator now has discretization errors of 

We note in passing that for this model one may introduce a topological 
charge in its 'geometric' definition z/ oc J2ti,uF^,u f^- Thus Ssw introduces a 
local coupling of to z/cra. 

The coefficient csw, which generally has to be determined non-perturba- 
tively at finite P, is just csw = 1 in 2D, since the gauge couphng e is di- 
mensionful and therefore csw = 1 + 0(e^ a^), which only gives corrections to 
higher order in a. 

Since the gauge part of the Wilson action 

Sg = ^T^MUp + uD (5) 

has C(a^) discretization errors, the total action 

S = Sw + Ssw + Sg (6) 

together with ^ gives a lattice version of the Schwinger model with dis- 
cretization errors of (9(a^). 
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We have done a MC simulation of the Schwinger model, using both 
the Wilson and the SW action, with two degenerate flavors of dynamical 
fermions. We used a hybrid Monte Carlo algorithm with 10 integration steps 
and a trajectory length tuned for a 0.8 acceptance rate of the MC step. For 
fermion matrix inversion we used a standard conjugate gradient algorithm, 
and there was no problem with its convergence or stability. In that algo- 
rithm the doubling of the fermionic species is required in order to guarantee 
positivity of the fermionic determinant - one of our reasons to study the 
Schwinger model with two flavors. 

In our implementation, the speed loss when including the clover term into 
the action is consistently ~ 30%. 

Depending on the lattice size (from 8 x 8 up to 18 x 18) and couplings, 
we skipped 20 to 100 configurations between measurements and finally had 
between 2 x 10'^ and 10^ independent configurations. To check for ergodicity, 
we measured the topological charge (using the geometric definition) and ob- 
served, that in each run the system tunneled several times between different 
topological sectors. 

3 Results 

3.1 Phase diagram 

The additional term in the action affects the critical n value where the theory 
has effectively massless fermions. Since the clover action is supposedly closer 
to the continuum (has smaller corrections to leading scaling behavior), we 
might expect less renormalization of with a value closer to its continuum 
value of 0.25 even for small 13. 

In order to determine Kc{(3) we follow suggestions [7^ to utilize the PCAC 
relation as discussed in ^8] for the Schwinger model. At each (3 we measured 
an observable proportional to the effective fermion mass for 5 different k 
values and determined Kc- Our result for a 8 x 8 lattice is compared with the 
result for the original Wilson fermion action in Fig. ^ 

The values k^. for the clover action are in fact closer to 0.25. However, 
both, Wilson and clover action have corrections 

= 0.25 + 0(0^) = 0.25 + C(l//3) . (7) 

The dependence on the finite lattice volume is weak and - for the Wilson 
action - is discussed elsewhere |H] . 
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Figure 1: Phase diagram for Wilson and clover action on an 8 x 8 lattice. 
The diamonds indicates the results (Wilson action) of an extrapolation to 
infinite lattice volume from results at various lattice sizes |H]. We find that 
Kf. is closer to its continuum value 0.25 for the clover action. 



3.2 Spectrum of the Dirac operator 

The spectral distribution of the Dirac operator is of some interest with re- 
gard to topological and chiral properties of the system. Recent studies for 
the Schwinger model [HI EH El Ej and for the SW-action in 4D ^3] have 
emphasized these aspects. Fig. |2] should be compared with typical spectra 
for the pure Wilson Dirac operator (cf. e.g. 

As compared to the Wilson action with its A A and (1 — A) ^ (A — 1) 
symmetries we observe for the SW operator an agglomeration of eigenvalues 
at larger distances from 0. This feature may improve somewhat the sepa- 
ration of the low-lying states from the doubling states already at moderate 
13. For the continuum limit the distribution density of small eigenvalues is of 
relevance. However, the distribution is not much closer to the circular shape 
one obtains for fixed point actions [Hj or Neuberger-projected actions [T5] . 
We cannot identify a clear signal of improvement in this respect for the SW 
Dirac operator. 
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Figure 2: The typical shape of the eigenvalue spectrum for the SW action; 
here shown for a configuration on a 16 x 16 lattice determined at /5 = 2 and 



3.3 Rotation invar iance 

We measured the propagators of the pseudoscalar isotriplet (vr) mesonic 
bound states. In the 2-fiavor model these modes are expected to be mass- 
less at Kf.. For short distances on the lattice, the rotation invariance of these 
propagators is clearly broken when using Wilson action. Generally, improved 
actions are expected to show better rotation invariance (cf. the Schwinger 
model FP action study in [16]). 

One must be careful not to compare these quantities simply at the same 
values of (3 and k, but at values corresponding to comparable discretization 
of the same continuum theory. We choose to compare the propagators at 
(5 = 2 and the respective Kc, so that we have a discretization of the massless 
(fermion) theory in both cases. 

The results are shown in FiglHl We find no noticeable improvement of 
rotation invariance when using the clover action. The relative error of the 
diagonal (1, 1) propagator compared to the (1, 0) and (2, 0) propagators is 
~ (18 ±1)% for both actions. This is the first indication, that short distance 
(high momentum) observables are not improved. 
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distance x in lattice spacings 

Figure 3: Plot of vr propagator as measured at /? = 2 on a 16 x 16 lattice 
at the respective k^. Using the clover action gives no improvement in the 
rotation invariance of the propagator. 

3.4 Meson dispersion relations 

Fig. 0] shows the dispersion relations for the massive (isosinglet vector) and 
the massless (isotriplet vector) meson at /3 = 2 and k = on a 12 x 12 
lattice. No improvement for the high momentum behavior is observed. 

Concerning the mass of the massive boson, however, the improvement is 
significant. Fig. El shows the ratio of the observed Schwinger mass over its 
continuum value at /3 = 2, k = on a 12x 12 and (3 = 4.5, k = on a 18x 18 
lattice. Assuming the continuum scaling behavior a(/3) oc these are 

different discretizations of the same physical system size La ex. Lj^f^ ~ 8.5. 
The improvement when including the clover term is quite obvious. The 
scaling corrections are remarkably smaller for the clover results. 

This behavior, the considerable improvement of low momentum states 
and the lack thereof for high momentum states, is no surprise when one con- 
siders, that the construction of Symanzik improved actions aims at improving 
the fermion dispersion relations for low momenta (cf. PP). The overall im- 
provement of the fixed point action in jTH] (which, however has more than 
hundred terms per site in the action) is definitely more pronounced and ex- 
tends to high momenta as well. 
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Figure 4: Dispersion relation (for = 2, lattice size 12 x 12) for a massless 
(isotriplet-vector) and the massive (isosinglet-vector) meson. The dashed 
curves denote the continuum dispersion relations. High momentum values 
are not improved by the clover action. 
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Figure 5: The Schwinger mass as measured at /3 = 2 on a 12 x 12 and 
/3 = 4.5 on a 18 X 18 lattice over the continuum Schwinger mass. For this low 
momentum observable, the improvement is significant. The lattice spacing 
is given in multiples of the value at /5 = 2, i.e. J'2./ j3. 
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4 Conclusion 



We performed a Monte Carlo simulation on the 2-flavor lattice Schwinger 
model, comparing standard Wilson to an 0{a) improved Sheikholeslami- 
Wohlert action. We found the following results: 

• The critical value of the hopping parameter Kc moves closer to its con- 
tinuum value 0.25. 

• The eigenvalue spectrum of the lattice SW Dirac operator changes but 
does not appear to improve significantly towards, e.g. a circular shape. 

• There was no improvement in the rotation invariance of meson propa- 
gators. For the meson dispersion relations, there is no improvement in 
the high momentum states. 

• There is considerable improvement in the scaling behavior of the Schwin- 
ger mass of the massive boson. 

These observations are consistent with the fact that the SW action is 
constructed as an improved action for low momentum fermionic states p <^ 
1/a. While there is considerable improvement in low momentum observables, 
high momentum observables are largely unaffected by the addition of the 
clover term. 

For the full fermion QCD simulations this suggests that it is preferable to 
use the clover action as long as one is only interested in the mass spectrum 
of the theory. On the other hand, some observables will hkely show no 
improvement at all, especially those connected with high momentum states. 
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